A comparison study of medium-modified QCD shower evolution scenarios 
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The computation of hard processes in hadronic collisions is a major success of perturbative Quan- 
tum Chromodynamics (pQCD). In such processes, pQCD not only predicts the hard reaction itself, 
but also the subsequent evolution in terms of parton branching and radiation, leading to a parton 
shower and ultimately to an observable jet of hadrons. If the hard process occurs in a heavy-ion 
collision, a large part of this evolution takes place in the soft medium created along with the hard 
reaction. An observation of jets in heavy-ion collision thus allows a study of medium-modified QCD 
shower evolution. In vacuum, Monte-Carlo (MC) simulations are well established tools to describe 
such showers. For jet studies in heavy-ion collisions, MC models for in-medium showers are currently 
being developed. However, the shower-medium interaction depends on the nature of the microscopic 
degrees of freedom of the medium created in a heavy-ion collision which is the very object one would 
like to investigate. This paper presents a study in comparison between three different possible im- 
plementations for the shower-medium interaction, two of them based on medium-induced pQCD 
radiation, one of them a medium-induced drag force, and shows for which observables differences 
between the three scenarios become visible. We find that while single hadron observables such as 
Raa are incapable of differentiating between the scenarios, jet observables such as the longiudinal 
momentum spectrum of hadrons in the jet show the potential to do so. 

PACS numbers: 25.75.-q,25.75.Gz 



I. INTRODUCTION 



Jet quenching, i.e. the energy loss of hard partons cre- 
ated in the first moments of a heavy ion collision due to 
interactions with the surrounding soft medium has long 
been regarded a promising tool to study properties of 
the soft medium H H, H, 0, [E Hi • The basic idea is to 
study the changes induced by the medium to a hard pro- 
cess which is well-known from p-p collisions. A number 
of observables is available for this purpose, among them 
suppression in single inclusive hard hadron spectra Raa 
i7:|, the suppression of back-to-back correlations [1, Q, 
single hadron suppression as a function of the emission 
angle with the reaction plane [1(3] and most recently also 
preliminary measurements of jets have become available 

Single hadron observables and back-to-back correlations 
are well described in detailed model calculations us- 
ing the concept of energy loss [H, 0, G3 , i-e- under 
the assumption that the process can be described by a 
medium-induced shift of the leading parton energy by 
an amount AE, followed by a fragmentation process us- 
ing vacuum fragmentation of a parton with the reduced 
energy. However, there are also calculations for these ob- 
servables in which the evolution of the whole in-medium 
parton shower is followed in an analytic way [H, [n| [ijj • 
Recently, also Monte Carlo (MC) codes for in-medium 
shower evolution have become available [H, EH US [HI 
which are based on MC shower simulations developed for 
hadronic collisions, such as PYTHIA [H or HERWIG 



'Electronic address: trenk@phys.jyu.fi 



[23|. In medium-modified shower computations, energy 
is not simply lost but redistributed in a characteristic 
way. 

All current in-medium shower MC codes model the in- 
teraction between partons and the medium in a different 
way. JEWEL (Jet Evolution With Energy Loss) [H| as- 
sumes either elastic collisions with thermal quasiparticles 
or, to implement radiative energy loss, an enhancement 
of the singular part of the parton branching kernels. Ya- 
JEM (Yet another Jet Energy- loss Model) [ll, i(| makes 
the assumption that the virtuality of partons traversing 
the medium grows according to the medium transport 
coefficient q which measures the virtuality gain per unit 
pathlcngth, and this medium-induced virtuality leads to 
increased radiation. Finally, Q-PYTHIA, the code pre- 
sented in [2l[ is a direct extension of the leading parton 
energy loss computations done in [1, [24[ and uses the dif- 
ferential radiation probabilities originally computed from 
a single hard parton now for each parton propagating in 
the shower simulation. At this stage, it is hardly sur- 
prising that different models employ different implemen- 
tations of the parton-medium interaction, as the nature 
of this interaction crucially depends on the microscopic 
properties of the medium, i.e. the very thing one wishes 
to determine from the experiments. 

A suitable strategy to determine these properties is thus 
to study the effects of various different implementations 
of the parton-medium interaction for different observ- 
ables. In this paper, we begin such a program by inves- 
tigating the effects on a number of different observables 
resulting from three different scenarios: Medium-induced 
radiation by an increase of parton virtuality dependent 
on the medium q as used in [Til, [2(j, an enhancement 
of the singular parts of the branching kernel leading to 
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additional radiation as used in [18|, [25[ and a drag force. 
Momentum-dependent drag forces appear in computa- 
tions modelling QCD-like N — 4 super Yang-Mills the- 
ories via the AdS/CFT conjecture [2g|, in the present 
paper we use a simplified ansatz in which a parton in 
a constant medium undergoes a momentum independent 
energy loss per unit pathlcngth. Such a drag term has 
not been tested in an in-medium shower evolution MC 
code previously. 

The paper is organized as follows: First, we briefly review 
the computation of medium-modified hadron jet as done 
in [l9|, [2(| • In addition, we describe the three different 
implementations of the parton-medium interaction and 
its relation to the spacetime structure of the shower in 
detail. In a first comparison, we make the connection to 
previous leading parton energy loss calculations by con- 
sidering a constant medium with fixed length. In this 
medium, we study the energy loss of the leading parton 
and present the result in terms of energy loss probabil- 
ity distributions and mean energy loss as a function of 
the parameters characterizing the medium. In a second 
comparison, we turn to a medium model which is closer 
to the experimental situation in so far as it expands and 
hence dilutes as a function of time. We compute various 
jet observables in this medium, such as the longitudi- 
nal momentum distribuion inside the jet or the angular 
broadening. Finally, in a last comparison we compute 
(as done in [l9j]) the suppression of the inclusive single 
hard hadron spectrum in terms of the nuclear suppression 
factor Raa and compare all scenarios with experimental 
results 0| • From this comparison, we tentatively deduce 
the relevant medium parameters. We conclude with a 
discussion of the implications of the results. 



II. MEDIUM-MODIFIED SHOWER 
EVOLUTION 



part of PYTHIA. 



A. Shower evolution in vacuum 

We model the evolution from some initial, highly vir- 
tual parton to a final state parton shower as a series of 
branching processes a — > b+c where a is called the parent 
parton and b and c arc referred to as daughters. In QCD, 
the allowed branching processes are q — ► qg, g — ► gg and 
g — > qq. The kinematics of a branching is described in 
terms of the virtuality scale Q 2 and of the energy fraction 
z, where the energy of daughter b is given by Ef, — zE a 
and of the daughter c by E c — (1 — z)E a . It is convenient 
to introduce t = InQ 2 / Aq C d where Aq C d is the scale 
parameter of QCD. t takes a role similar to a time in the 
evolution equations, as it describes the evolution from 
some high initial virtuality Qo (to) to a lower virtuality 
Qm (t-m) at which the next branching occurs. In terms 
of the two variables, the differential probability dP a for 
a parton a to branch is [2^ . |3C| 



dPa = Y l ^ L Pa^bc(z)dtdz 
' — ' Z7T 



b,c 



(1) 



where a s is the strong coupling and the splitting kernels 
P a -> bc (z) read 



P q ^ qg (z)=A/3^ 



Pg^gg( z ) — 3 



(l-z(l~z)f 



(2) 
(3) 



P g ^ qq (z) = N F /2(z 2 + (f - zf) (4) 



In this section, we describe how the medium-modified 
fragmentation function (MMFF) is obtained from a com- 
putation of an in-medium shower followed by hadroniza- 
tion. Key ingredient for this computation is a pQCD MC 
shower algorithm. In this work, we employ a modifica- 
tion of the PYTHIA shower algorithm P YSHOW . In 
the absence of any medium effects, our algorithm there- 
fore corresponds directly to the PYTHIA shower. Fur- 
thermore, the subsequent hadronization of the shower is 
assumed to take place outside of the medium, even if the 
shower itself was medium-modified. It is computed using 
the Lund string fragmentation scheme [28| which is also 



where we do not consider electromagnetic branchings. 
Np counts the number of active quark flavours for given 
virtuality. 

At a given value of the scale t, the differential probability 
for a branching to occur is given by the integral over all 
allowed values of z in the branching kernel as 



[ z+(t) a s 

Ia->bc(t) = / dz-^P a ^ bc (z). 
Jz-(t) 2lr 



(5) 



The kincmatically allowed range of z is given by 
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4MJMI 



(6) 



where Mf = Q 2 +m 2 with n%i the bare quark mass or zero in the case of a gluon. Given the initial parent virtuality 
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Q\ or equivalently t a , the virtuality at which the next 
branching occurs can be determined with the help of the 
Sudakov form factor S a (t), i.e. the probability that no 
branching occurs between to and t m , where 



S a (t) = exp 



b,t 



(7) 



Thus, the probability density that a branching of a occurs 
at t m is given by 



dPg 

dt 



b.c 



S a (t). 



(8) 



These equations are solved for each branching by the 
PYSHOW algorithm [27| iteratively to generate a shower. 
For each branching first Eq. {§} is solved to determine the 
scale of the next branching, then Eqs. ©-© are evalu- 
ated to determine the type of branching and the value of 
z, if the value of z is outside the kinematic bound given 
by Eq. |6]) then the event is rejected. Given to,t m and 
z, energy- momentum conservation determines the rest 
of the kinematics except for a radial angle by which the 
plane spanned by the vectors of the daughter parents can 
be rotated. 

In order to account in a schematic way for higher order 
interference terms, angular ordering is enforced onto the 
shower, i.e. opening angles spanned between daughter 
pairs &, c from a parent a are enforced to decrease ac- 
cording to the condition 



z b (l - z)b) 



Mi 



> 



(9) 



After a branching process has been computed, the same 
algorithm is applied to the two daughter partons treating 
them as new mothers. The branching is continued down 
to a scale Q m in which is set to 1 GeV in the MC simula- 
tion, after which the partons are set on-shell, adjusting 
transverse momentum to ensure energy-momentum con- 
servation. 

After all possible branchings have been performed, i.e. 
after for all partons the condition Q < Q m in has been 
reached, the resulting parton shower is connected with 
a string following the Lund scheme [28| which is subse- 
quently allowed to decay into hadrons. These hadrons 
form the observable jet, and analyzing the distribution 
of hadrons, we may for example determine the fragmen- 
tation function Df^h{z), i.e. the distribution of hadron 
species h with an energy Eh — zEf originating from a 
shower initiating parton / where Ef is the whole energy 
of the jet. 



B. Spacetime structure of the shower 

While the vacuum shower evolution equations above are 
solved in momentum space only, the interaction with the 
medium requires modelling of the shower evolution in 
position space as well, because the medium properties 
in a general medium change as a function of the posi- 
tion space variables. Usually, these are given in the cm. 
frame of the collision in terms of the spacetime rapidity 
r) s , the radius r, the proper time r and the angle </>, and 
knowledge of the medium evolution implies knowledge of 
medium properties such as the local medium temperature 
T in the form T(j] s , r, 0, t). 

In order to make the link from momentum space to mo- 
mentum space, we assume that the average formation 
time of a shower parton with virtuality Q is developed 
on the timescale 1/Q, i.e. the average lifetime of a virtual 
parton with virtuality Qb coming from a parent parton 
with virtuality Q a is in the rest frame of the original hard 
collision (the local rest frame of the medium may be dif- 
ferent by a flow boost as the medium may not be static) 
given by 



(n) 



Ql 



Eb 
r>2 



(10) 



Going beyond the ansatz of [19|, |20( where we used this 
average formation time for all partons, in the present 
work we assume that the actual formation time can be 
obtained from a probability distribution 



P(n) = exp 



n 
(n) 



(li) 



which we sample to determine the actual formation time 
of the fluctuation in each branching. This establishes 
the temporal structure of the shower. With regard to 
the spatial structure, we make the simplifying assump- 
tion that all partons probe the medium along the cikonal 
trajectory of the shower initiating parton, i.e. we neglect 
the small difference of the velocity of massive partons to 
the speed of light and possible (equally small) changes 
of medium properties within the spread of the shower 
partons transverse to its axis. 



C. The parton- medium interaction 

In the following, we assume that any effect of the medium 
will affect the partonic stage of the evolution, but not 
the hadronization. This is equivalent to the idea that 
hadronization takes place outside the medium, an as- 
sumption commonly made also for leading parton energy 
loss calculations. The validity of this assumption will be 
dicussed below. 

We use three different scenarios to model the interaction 
of partons with the medium. The first one, in the fol- 
lowing referred to as RAD, has been used previously in 



4 



[T^ . [20j . The relevant property of the medium probed 
is the transport coefficient q(r) s ,r,(f),T) which represents 
the virtuality gain AQ 2 per unit pathlength of a par- 
ton traversing the medium. Note that this represents 
an average transfer, i.e. a picture which would be re- 
alized in a medium which is characterized by multiple 
soft scatterings with the hard parton. However, unlike 
in p^ . [2tjj the virtuality transfer to a shower parton is 
randomized in the present work since the formation time 
is distributed randomly around its average. Thus, effec- 
tively the present scenario includes the possibility to have 
both a small formation time and hence a small virtuality 
gain and a large formation time corresponding to a more 
substantial increase in virtuality. 

In practice, we increase the virtuality of a shower par- 
ton a propagating through a medium with specified 
q(Vs,r,<j>,T) by 

AQ 2 , = / d(q(C) (12) 

where the time r a is given by Eq. (|lip . the time t° is 
known in the simulation as the endpoint of the previous 
branching process and the integration d£ is along the 
cikonal trajectory of the shower-initiating parton. If the 
parton is a gluon, the virtuality transfer from the medium 
is increased by the ratio of their Casimir color factors, 
3/| = 2.25. 

If AQ 2 -c Q 2 , holds, i.e. the virtuality picked up from 
the medium is a correction to the initial parton virtuality, 
we may add AQ 2 to the virtuality of parton a before 
using Eq. ([8]) to determine the kinematics of the next 
branching. If the condition is not fulfilled, the lifetime 
is determined by Q 2 + AQ 2 and may be significantly 
shortened by virtuality picked up from the medium. In 
this case we iterate Eqs. (ITU)) . (fT2|) to determine a self- 
consistent pair of ((r a ), AQ 2 ). This ensures that on the 
level of averages, the lifetime is treated consistently with 
the virtuality picked up from the medium. The actual 
lifetime is still determined by Eq. (fTTj) . 
In a second scenario, in the following called DRAG, we 
assume that the medium exerts a drag force on each prop- 
agating parton. The medium is thus characterized by a 
drag coefficient D(r] s ,r,4>,T) which describes the energy 
loss per unit pathlength. 

In the simulation, the energy (and momentum) are re- 
duced by 

f r°+r a 

AE a = / d(D(() (13) 
J T ° 

Again, for a gluon the energy loss is increased by the color 
factor ratio 2.25. As in the previous case, the energy loss 
induced by the drag force is randomized even given the 
branching kinematics due to the randomized formation 
time of a branching. 



The third scenario has been suggested in [l8|, [2|| . In the 
following, it is referred to as FMED. Here, the modifica- 
tion does not concern the parton kinematics, but rather 
the evolution kernel, Eqs. (J2H11)- ln this scenario, the 
singular part of the branching kernel in the medium is 
enhanced by a factor 1 + f me d, e.g. Eq. © becomes in 
the medium 

(14) 

The effect of the medium is thus summarized in the value 
of fmed- Note that in the FMED scenario, no explicit ref- 
erence to the spacetime structure of the shower is made, 
in this sense, the scenario is rather different from the 
other two. 

Note that in the RAD scenario the shower gains energy 
from the medium by means of the virtuality increase, 
in the DRAG scenario the shower loses energy to the 
medium whereas the shower energy is conserved in the 
FMED scenario. While this appears surprising at first, 
it is actually rather a matter of book-keeping. For a 
shower in the medium, there is no conceptual way to sep- 
arate soft partons from the shower and from the medium. 
However, the model framework outlined above does not 
treat the medium as consisting of partons, but rather as 
an effective influence on the shower. Thus, in a more 
realistic model one would define a criterion (say a mo- 
mentum scale) based on which partons are removed from 
the shower and become part of the medium. In such a 
model, all three scenarios would lead to a loss of energy 
from the shower to the medium through the appearance 
of soft partons in the evolution, in addition to possible 
other mechanisms of energy transfer to the medium. 



III. COMPARISON FOR A CONSTANT 
MEDIUM 

In this section, we perform several computations for the 
simple case of a constant medium with fixed pathlength. 
This is chiefly done in order to establish the relation of 
the models outlined above to older computations based 
on leading parton energy loss. 



A. Presence and absence of scaling 

A constant medium corresponds to a choice of a single 
value of g, D or f me d- However, in the case of both the 
RAD and the DRAG scenario, also the medium length L 
has to be specified, thus in principle the medium is char- 
acterized by two parameters. In [l!| however we found 
an approximate scaling law for the RAD scenario accord- 
ing to which the modification chiefly depends on the vir- 
tuality picked up along the eikonal path of the shower 
initiating parton AQ 2 ot = j d(q(() or in the case of a 
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FIG. 1: The MMFF of a d-quark into charged hadrons for constant value of AQ 2 ot — qL = 5 GeV 2 in the RAD scenario (left 
panel) and AEtot = DL = 5 GeV in the DRAG scenario (right panel) for different pathlengths in a constant medium. 



constant medium simply qL. A similar scaling law can 
also be established for the DRAG scenario, albeit only in 
the case of an expanding medium (see below) . Whenever 
such a scaling law holds, a comparison between the differ- 
ent scenarios can be made based on the single parameter 
AQtot or AE tot only. 

It is clear that the scaling cannot work for all the pos- 
sible functional forms q((). In two different limits this 
can be made plausible: If, in the RAD scenario, AQ\ ot is 
added at once initially, AQ 2 /Q 2 is for reasonable values 
of hard process kinematics and medium properties very 
small. For example, for typical RHIC kinematics the ini- 
tial Q 2 from which the evolution starts may be 400 GeV 2 
whereas the total virtuality acquired for a parton travers- 
ing the whole medium is about 15 GeV 2 according to the 
results of [lj|. However, such a small correction will not 
influence the shower evolution significantly. On the other 
hand, if the virtuality is added later when the typical Q 2 
is of order of AQ 2 , a much stronger modification is ex- 
pected. Thus, one expects the scaling law to be violated 
into the direction of less medium effect if q(Q is strongly 
peaked towards r = 0. 

A similar argument, can be made for the DRAG scenario. 
The drag force acts on every parton in the shower. This 
means that if D(() is strongly peaked towards r = 0, 
then the drag force acts only on one parton, the shower 
initiator, whereas if it is applied later, its effect is felt by 
several partons. 

On the other hand, note that the shower evolution is 
terminated for every parton which reaches Q 2 < Q 2 ^ = 
1 GeV 2 . This implies that the typical lifetime of the 
shower for an initial parton with energy E is given by 
Tmax ~ EjQ 2 minl thus a shower with E = 20 GeV probes 
the medium on average for a distance of 4 fm (Eq. (fTTj) 
leads to fluctuations around this average though). Thus, 
if L is chosen much beyond T max , qL or DL are not good 
parameters any more, as the shower does not effectively 
probe the whole medium. 



The latter effect is clearly not related to an actual physics 
effect but rather an artefact of the need to switch to 
a non-perturbative description of hadronization at some 
point in the simulation. It is unreasonable that a parton 
(or proto-hadron) would feel no effect from the medium 
just because its virtuality is small, however it is unclear 
just how the effect should be implemented properly in 
the present framework. The behaviour of the simulation 
thus depends on the actual choice of Q m in, and this needs 
to be optimized eventually in comparison with data. A 
study of the effect of changing Q mm will be presented 
below. 

The resulting MMFF for a light quark into charged 
hadrons for constant qL or DL and a variation of path- 
length is shown in Fig. [T] for both the RAD and the 
DRAG scenario. In a constant medium, the RAD sce- 
nario shows approximate scaling for pathlcngth between 
0.5 and 5 fm. The DRAG scenario does not exhibit a 
strong scaling in the region of large z, but in the region 
z ~ 0.5 which is predominantly probed when computing 
the single hadron spectra, the variations are not too large 
for pathlengths between 0.5 and 3 fm. 

Note that neither the short pathlength nor the long path- 
lcngth limit is actually problematic for a realistic medium 
evolution taken from a hydrodynamical model. The first 
limit is avoided by virtue of the thermalization time of or- 
der O(0.6) fm for RHIC kinematics. This is a large time 
compared with the timescale in which the first branch- 
ings in a shower occur, thus by the time the medium is 
present, the shower is already well developed. The sec- 
ond limit is avoided because in an expanding medium 
q{Q or D(() drop rapidly as a function of time, thus the 
late time contribution to J e2C?(C) or / ^C-D(C) is small 
in any case. Thus, the scaling works much better for a 
realistic evolution as the constant medium results would 
suggest. 
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B. Energy loss and quenching weights 

We now proceed to compare the three scenarios on the 
basis of leading parton energy loss. This is relevant 
to make the connection to previous calculations in the 
BDMPS or ASW formalism [2, HJ] which are formulated 
using this concept. For this purpose, we select the shower 
initator to be a c-quark and extract the energy distri- 
bution of the leading c-quark dN/dE c after the shower. 
From the comparison of the distribution dN/dE^ ac in 
vacuum and in the medium dN/dE™ ed , we can deduce 
the energy loss probability distribution P(AE). The idea 
is to make an ansatz 



is very different for a light quark shower where the lead- 
ing quark distribution typically peaks at z ~ 0.5 and any 
energy loss of AE > E/2 shifts the bulk of the distribu- 
tion into the unphysical region of negative energies. 

In Fig. [2] we show the leading charm distributions both in 
vacuum and in medium for a medium pathlength of L = 2 
fm. In order to make a meaningful comparison between 
the different scenarios, the average relative energy loss 
(AE)/E is fixed to 10% or 20% respectively. 

In order to deduce the energy loss probability distribu- 
tion from these results, we have to solve Eq. (fT5|) . By 
discretizing the integral over AE in Eq. (p~5|) we can cast 
it into the form of a matrix equation 



dN 
dE, 



med n j at- vac 

(E) = / d{AE) — {E')P(AE)S{E' - E - AE) 
J dE c 

(15) 

and solve it for P{AE). Note that this ansatz contains 
the rather drastic assumption that there is no parametric 
dependence on the initial energy E. If we require P(AE) 
to be a probability distribution, the assumption may im- 
ply that for some partons in the distribution dNjdE v c ac 
the energy loss AE is larger than their energy E in which 
case they have to be considered lost to the medium. A 
similar situation also occurs in the application of the 
ASW formalism to finite energy kinematics. The problem 
of the validity of assuming energy independence however 
only concern the comparison with the ASW results in 
which energy loss is formulated in terms of a probabil- 
ity density P(AE). In all other results presented in this 
manuscript, the full information of the shower including 
finite energy kinematics is used and no assumption about 
energy independence of energy loss needs to be made. 
The choice of a c-quark as shower initiator has a twofold 
motivation. First, it allows to define energy loss in the 
same way as done in the ASW formalism. Note that the 
ASW formalism assumes infinite parent parton energy 
and calculates energy loss via the radiation spectrum off 
the parent. In applying the formalism to finite energy, a 
process may occur in which a radiated gluon takes 90% 
of the energy of an initial quark q\ . This energy is then 
considered to be lost from the q±. However, in the shower 
language, the radiated gluon would in this case become 
the new leading parton, and even tagging the leading 
quark out of a shower would not prevent processes where 
this gluon splits into a qq pair where the new quark q2 
might still be harder as the original parent q\ of the gluon. 
The choice of a c quark as shower initiator effectively 
suppresses such processes and allows to treat energy loss 
as closely as possible to ASW [31 1. 

The second advantage of extracting P(AE) from a c 
quark is that the c-fragmentation is rather hard, i.e. the 
probability distribution to find the leading c-quark af- 
ter a vacuum shower peaks close to 2 = 1. This effec- 
tively means that if one considers an additional, medium- 
induced shift in energy, most of the energy range is still 
available for the dominant part of the distribution. This 



Ni(E i ) =Y J K ij {E\AEi)P j (A&) 



(16) 



where dN/dE c is provided at to discrete values of E la- 
belled Ni and P(AE) is probed at n discrete values of 
AE labelled Pj. The kernel Kij is then the calculated 
dN/dE c for all pairs (E l ,AEi) where the energy loss 
acts as a shift of the distribution, i.e. dN/dE™ ed (E) — 
dN/dEl ac (E + AE). 

Eq. (TT6")) can in principle be solved for the vector Pj by 
inversion of for m = n. However, in general this does 
not guarantee that the result is a probability distribution. 
Especially in the face of statistical errors and finite nu- 
merical accuracy the direct matrix inversion may permit 
negative Pj which have no probabilistic interpretation. 

Thus, a more promising solution which avoids the above 
problems is to let to > n and find the vector P which min- 
imizes \\N — KP\\ 2 subject to the constraints < Pi < 1 
and X)"=i -Pi = !■ This guarantees that the outcome can 
be interpreted as a probability distribution and since the 
system of equations is overdetermined for m > n errors 
on individual points Ri do not have a critical influence 
on the outcome any more. This is the approach we have 
chosen. 

The results are shown for L = 2 fin in Fig. [3J Qual- 
itatively, both the radiative energyloss scenarios RAD 
and FMED produce energy loss probability distributions 
which are similar to the ASW quenching weights [24| in 
the sense that they are flat across a wide range in AE. In 
contrast, the DRAG scenario produces a localized peak in 
the energy loss distribution which reminds of the quench- 
ing weights found for elastic energy loss scenarios [32l[33|. 
Especially for larger energy loss, the RAD and the FMED 
scenario lead to almost identical results. 

However, there is an important difference to the ASW 
quenching weights: While the ASW results typically 
show a large discrete probability for no energy loss, the 
results obtained here show no substantial strength in the 
first bin (the inversion procedure outlined above cannot 
separate zero energy loss from small energy loss). 
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FIG. 2: Energy distribution of the leading c quark for a 20 GeV c quark as shower initiator in the three different scenarios for 
the parton-medium interaction (see text). Left panel: 10% average energy loss, right panel: 20% average energy loss. 
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FIG. 3: Energy loss probability distribution P(AE) for the leading c-quark for a 20 GeV c-quark as shower initiator in the 
three different scenarios for the parton-medium interaction (see text). Left panel: 10% average energy loss, right panel: 20% 
average energy loss. 



C. Parametric dependence of mean energy loss 

In order to gain more insight into the different sce- 
narios, we investigate in Fig. U for a constant medium 
with L = 2 fm how the mean energy loss, defined as 
(AE) = J dAEAEP{AE) with P(AE) obtained as in 
the previous section behaves as a function of the relevant 
medium parameters. We include a scenario in which the 
strong coupling constant is not allowed to run with the 
virtuality scale in the shower (as is the default option in 
PYSHOW) but is kept fixed at a s = 0.3. 
There is no unique way to present and compare the re- 
sults, as the three relevant parameters q, D and fmed are 
rather different. However, as apparent from Fig. [H it is 
possible to find a simple proportionality relation between 
q and f me d such that the rise of the mean energy loss ap- 
pears very similar. This, in addition to the similarity of 
P(AE) for both the RAD and the FMED scenario points 



towards some generic properties of radiative energy loss 
scenarios independent of the details of the implementa- 
tion. 



In particular, the RAD and the FMED scenario exhibit 
saturation of the mean induced energy loss at about 25% 
of the total energy as the medium effect is increased. 
This saturation is even more pronounced for a constent 
a s . In striking contrast, the DRAG scenario in which 
energy is directly transferred to the medium shows an 
almost linear rise up to mean energy losses of 50%. Note 
that the extraction of the energy loss probability based 
on discretization and matrix inversion as outlined above 
becomes increasingly problematic at (AE) /E > 0.4 due 
to the problem of partons being shifted to negative ener- 
gies mentioned above. 
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FIG. 4: Mean energy loss as a function of the medium prop- 
erties in different scenarios for the parton-medium interaction 
in a constant medium with L = 2 fm. 



IV. COMPARISON FOR A SINGLE PATH IN 
AN EXPANDING MEDIUM 

We now turn to a more realistic scenario in which the 
parton propagates in a medium as created in a heavy- 
ion collision. Rclativistic fluid-dynamical models such 
as [34| give a good description of many bulk proper- 
ties of the medium, hence in the following we will as- 
sume that hydrodynamics is a valid description of the 
medium. Both the finite size and the finite lifetime of 
such a medium are felt by the parton. In particular, the 
local density may drop a) because of a spatial variation, 
i.e. the parton reaches the medium edge and b) a tem- 
poral variation, i.e. the global expansion of the medium 
reduces the overall density as a function of time. In addi- 
tion, there are arguments that the hydrodynamical flow 
of the medium expansion should also have a direct influ- 
ence on the medium properties as seen by the medium 
due to Lorentz transformation between the moving local 
medium rest frame and the frame of the hard collision 




A. Characterization of the medium 

In we have established that if q is linked with the 
medium properties by the relation 

q(() = K-2- [e(C)] 3/4 (coshp(C) - sinhp(C) cost/;) (17) 

with K a parameter determining the interaction strength 
which is a priori unknown (in an ideal QGP, K — 1 is 
expected [35[ but a comparison study of different energy 
loss models has shown to be inconclusive in extracting 
values for K [36|). the medium energy density e, the lo- 
cal flow rapidity p with angle ip between flow and par- 



ton trajectory [37|, |38|, we find that the vast majority of 
paths found in the 3-dimcnsional hydrodynamical model 
of Bass and Nonaka [11] leads to a q(() which can be 
described by the rather simple expression 

* (0 = (b + r/dfm/c)y (18) 

Based on this expression, we investigated three different 
scenarios (approximately representing a parton travelling 
into +x direction originating from x = 4 fm (A), x = 

(B) and x — —4 fm (C), y = in all cases in the trans- 
verse (x, y) plane at midrapidity. These trajectories are 
characterized by the parameters (b — 1.5, c = 3.3, te = 
5.8 fm/c) (A), (b = 1.5, c = 2.2, t e = 10 fm/c) (B) and 
(b = 1.5, c = 2.2, te = 15 fm/c) (C) and are quite typ- 
ical for partons close to the surface (A), emerging from 
the central region (B) or traversing the whole medium 

(C) . As in IIII Al in the present paper for a constant 
medium, we found that an approximate scaling in which 
the medium effects did not depend on details of the tra- 
jectories (A), (B), or (C) but only on AQ^ ot = Jd(q((). 
The virtue of this scaling law is twofold: First, it allows to 
present the medium modifications for the relevant class of 
functions g(C) as a function of a single parameter AQt ot 
only. Second, it considerably speeds up the computation 
for a comparison with data where a weighted average 
over all possible paths through the medium has to be 
computed. 

In [19( we have made the rather drastic assumption that 
the medium does not exert any effect before the thermal- 
ization of the medium at the time Tj ra where Ti n = 0.6 
fm/c in the model studied for RHIC (34|. In the fol- 
lowing, we adopt a more realistic approach in which we 
increase the medium effect linearly from zero at r = 
to its value reached at Ti„. The idea behind this is that 
initially no medium can be present, as the timescale for 
hard processes precedes any other timescale in the sys- 
tem. However, even a medium which is not yet equili- 
brated may interact with hard partons and lead to scat- 
tering processes. A linear interpolation between the ini- 
tial time and the equilibration time seems a reasonable 
prescription to capture part of these effects. In practice, 
qualitative aspects of the results of [l9| . in particularly 
the presence of the scaling, are not substantially altered 
by this modification. There is however an effect on the 
numerical value of extracted medium parameters. 
Let us now consider the other scenarios DRAG and 
FMED. Eq. (fTTj) which links q with the hydrodynami- 
cal properties of the medium is based on counting the 
potential scattering centers along the parton trajectory. 
£ 3/4 f or an idggj g as corresponds to the entropy density, 
which in turn is proportional to the medium density. The 
additional factor (coshp(^) — sinh p(Q cosip) is nothing 
but the appropriate transformation to determine how the 
density seen by the parton is changed under a boost of 
the restframe of the medium [37j . It is reasonable to as- 
sume a similar measure of potential scattering centers to 
be relevant for the other scenarios. This ansatz leads to 
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D(f) = K D ■ [e(C)] 3/4 (coshp(C) - sinhp(C)cos^) (19) 

for the drag coefficient D with an a priori unknown pa- 
rameter Kjj specifying the overall strength of the drag 
force. 

As discussed above, the FMED scenario has no explicit 
dependence on the spacetime evolution of the shower, 
but it seems reasonable the the parameter f me d should 
depend on the total effect of the medium measured in the 
number of potential scatterers which have been encoun- 
tered. This leads to the ansatz 



f med = K f / rfC[e(C)] 3/4 (coshp(C) - smhp(C)cosV). 



(20) 

Here, as in the previous scenarios, we also introduce an 
a priori unknown parameter Kf which determines the 
strength of the parton-medium interaction. 
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FIG. 5: The MMFF of a 20 GeV d-quark into charged hadrons 
for three different paths (A), (B) and (C) (see text) with 
AEtot = 5 GeV in the DRAG scenario. 

The scaling within the RAD scenario of the results with 
AQ\ ot has been established in [l9[ and in the present 
paper also for a constant medium in IIII Al The DRAG 
scenario shows no strong scaling for a constant medium, 
but as anticipated the result is more promising for an 
expanding medium. The validity of the scaling under 
these conditions is apprent from Fig.[5]where we compute 
for fixed AE for the three different paths (A), (B) and 
(C). Note that scaling of the FMED scenario is realized 
by definition using the ansatz Eq. (f2"0"|) . 



B. Longitudinal momentum distribution of the 
shower 



In Fig. [5J we show the longitudinal momentum distri- 
bution of charged hadrons inside the shower in therms 



of the MMFF D(z) for three different scenarios in com- 
parison. There is no a priori criterion at which values 
of the three medium parameters AQf ot , AE to t and f me d 
the three different scenarios should be compared. For 
the comparison in terms of energy loss probability distri- 
butions done above we required a fixed value (AE)/E, 
but this is not a meaningful variable when one wants to 
compare on the basis of the whole parton shower instead 
of the leading parton kinematics only. Here, we chose 
the criterion that the MMFF approximately agree in an 
interval of 0.4 < z < 0.7. This is the region of the frag- 
mentation function which is predominantly probed when 
the fragmentation function is folded with a pQCD par- 
ton spectrum to compute single inclusive hadron produc- 
tion. The implication is that a computation with MMFFs 
agreeing in the above interval would yield approximately 
the same observable hadron spectra. This choice leads to 
the interesting and amusing numerical coincidence that if 
the parameters are given in powers of GeV, the relation 
AQ 2 tot /GeV 2 » A£/GeV « 10/ med holds. 
We show the MMFF of a 20 GeV d-quark into charged 
hadrons for AQ 2 = 10 GeV 2 (the parameters of the other 
scenarios adjusted correspondingly) in Fig.[6l right panel. 
In order to focus more on the hadron production at low 
momenta, we introduce the variable £ = ln(l/x) where 
x = p/Ej e t is the fraction of the jet momentum carried 
by a particular hadron and Ej et is the total energy of 
the jet. The inclusive distribution dN/d£, the so-called 
Hump-backed plateau, is an important feature of QCD 
radiation [39l |40] | and is in vacuum dominated by color 
coherence physics. 

In Fig. [5] (left panel) we show dN/d^ for the three differ- 
ent scenarios in comparison with the unmodified result. 
It is apparent from the figure that while the three scenar- 
ios agree in the high z and consequently low £ region, they 
exhibit sizeable differences in the high £ region where in- 
duced radiation is expected to contribute to soft hadron 
production. Here, both the radiative scenarios RAD and 
FMED show the expected enhancement of the distribu- 
tion, but the DRAG scenario is strikingly different — it 
falls below the vacuum result. However, this is hardly 
surprising, as in this scenario energy is taken away from 
the evolving shower and is hence not available for hadron 
production. 

While a measurement of dN/d^ would appear to be a 
promising means to distinguish between induced radia- 
tion and a drag force as the microscopic realization of 
energy loss, it has to be pointed out that there are two 
things which urge some caution. First, the Lund scheme 
used to model hadronization in the present framework 
assumes that hadronization takes place far outside the 
medium. If the energy of a hadron h of mass nih is Eh, 
the spatial scale at which hadronization occurs can be es- 
timated as lh ~ Eh/m^. For pions, this is not a problem 
throughout the kinematic range, but for kaons and pro- 
tons the hadronization length is considerably shortened. 
Even a 10 GeV proton has only lh ~ 2 fm, thus heavy 
hadron production in the high £ region is not addressed 
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FIG. 6: Longitudinal momentum distribution of charged hadrons inside a jet originating from a 20 GeV d-quark shown as 
fragmentation function D(z) (left panel) and dN/d£ (right panel) for the vacuum and the three different scenarios for the 
parton-medium interaction (see text). The medium parameters have been chosen to let the modified D(z) approximately agree 
for 0.4 < z < 0.7 for all in-medium scenarios. 



adequately in the model, as one cannot safely assume 
hadronization takes place outside the medium where the 
Lund model is applicable. Nevertheless, since pions con- 
stitute the bulk of charged hadron production, the essen- 
tial features of the model are expected to be robust. 
The second issue concerns the effect of trigger bias. A sc- 
ries of experimental cuts has to be imposed on events in 
heavy-ion collisions to discriminate hadrons belonging to 
jets from the background of soft medium hadrons. How- 
ever, strongly modified jets (for example those emerging 
from the medium center) are less likely to fall within the 
cuts than unmodified jets (such as those from the medium 
edge). As a result there is a trigger bias which suppresses 
events in which a modification of dN/dl; is visible. A cal- 
culation in the RAD scenario taking into account a re- 
alistic series of experimental cuts has been performed in 
[20] | and found that there should be no visible enhance- 
ment if jets are identified directly via a standard set of 
cuts. 



C. Angular distribution 

Another possibility to identify the mechanism of the 
parton-medium interaction is to study the structure of 
the jet transverse to the jet axis. This is reflected e.g. in 
the angular distribution of hadrons around the jet axis. 
The distribution dN/d<f> where 4> is t ne angle between 
hadron and jet axis for AQ 2 ot — 10 GeV 2 (the parameters 
in the other scenarios adjusted accordingly) for the vac- 
uum and the three different scenarios is shown in Fig. [7| 
where a cut in momentum of 1 GeV has been applied 
to focus on hadrons which would appear above the soft 
background of a heavy-ion collision. 
It is apparent from the figure that the radiative energy 
loss scenarios again roughly agree with each other and 



E = 20 GeV 

d -> h + ", P T > 1 GeV 




<|> [rad] 

FIG. 7: Angular distribution of charged hadrons above 1 GeV 
coming from the fragmentation of a 20 GeV d-quark for vac- 
uum and three different scenarios of parton-medium interac- 
tion (see text). 

lead to angular broadening of the jet as compared to 
the vacuum result, whereas the DRAG scenario shows 
no indication for broadening. 

D. The sensitivity to Q m in 

For a constant medium, we noted earlier that there is a 
sensitivity to the choice of the minimum virtuality scale 
Qmin at which partons in the shower are evolved further. 
Before comparing the results of this section to data, it 
is reasonable to ask to what extent a choice of Q m in 
different from its default value Q m in — 1 GeV has an 
influence on the results. 
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FIG. 8: MMFF for a 20 GeV d-quark as shower initiator for 
different values of the minimum shower virtuality Qmin in 
the RAD and the DRAG scenario (see text) where medium 
parameters have been adjusted to compensate for the choice 

of Q ■mm ■ 

In Fig. [U we show results for the MMFF in both the 
RAD and the DRAG scenario with a lower Qmin = 0.7 
GeV where the shower evolves on average a factor two 
longer. Superimposed are results with the default choice 
Qmin — 1 GeV for which the medium parameters AQ\ ot 
and AE to t respectively have been increased for the best 
possible agreement of the results. 

It is evident from the figure that a lower Q m in does not 
substantially influence the shape of the resulting MMFF, 
but that at least for RHIC kinematics, a lower choice of 
Qmin can be compensated by assuming a different choice 
of the medium parameters. It is thus not possible to 
extract definite values for q or D from a mode fit to mea- 
sured single hadron spectra, rather only pairs (q, Q m i n ) 
can be determined. 



V. COMPARISON WITH NUCLEAR 
SUPPRESSION DATA 

In this section, we aim at comparing with experimental 
data. This implies that neither initial position nor initial 
momentum nor type of the shower initiating parton are 
known. The probabilities to find a given parton type 
with given momentum have to be computed in pQCD 
whereas the probability to produce a parton at a given 
vertex position can be found from overlap calculations. 
Note that the need to average over position and initial 
momentum corresponds to a substantial increase in MC 
computing time which could not be done without using 
the scaling laws. 

A. The averaging procedure 

We begin the analysis by showing how to compute the nu- 
clear suppression factor R AA using the medium-modified 
fragmentation function in the hydrodynamically evolv- 
ing medium. For this, we first have to obtain the single 
inclusive hard hadron spectrum. 

We treat the partonic subprocesses of the hard reaction 
in leading order pQCD. The straightforward calculation 
involves the convolution of the initial nucleon [4l|, |42j (or 
nuclear [4^, 0, |45[) parton distribution functions with 
the relevant pQCD subprocesses and yields the single in- 
clusive distribution d<T , $ , of hard partons / in trans- 

verse momentum and rapidity yf where the rest of the 
reaction X is unobserved (more detailed expressions can 
be found in 0). 

The single inclusive hadron distribution in hadronic mo- 
mentum Pt and rapidity y follows from the parton spec- 
trum through the convolution with the fragmentation 
function D f^h(z, fi^) where z is the momentum fraction 
taken by the hadron and /i/ is the hadronic momentum 
scale as 



da AB^h+x r da AB ^f+ x f 1 



dP 2 dy 2^ J d PT d Vf dy\dy s J dzD s ^ h {z^))8{m 2 T - M^{ PT ,y f ,z))S(y~Y{p Tl y fl z)) (21) 



r 



with The nuclear suppression factor is defined as 

M*(PT,y f ,z) = { z PT f + MH^ yf , (22) R AA (P T ,y) = T^^kdy ^ 

2mT where T AA (h) is the standard nuclear overlap function. 

Zmin = ~/T cos ^ y We can compute it by forming the ratio 



and 

m T U1 J' K ' "*a\*i,i,j dp 2 dy I dp 2 dy 



P \ da AA ^ h+x d(T pp ^ h+x 

Vb T ,y / ,z)=arsinh(-^sinh J// . (24) R AA (P T ,y) = "i*"" / ( 26 ) 
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where da^ h + x jdP T dy follows from Eq. flUJ when 
Df^h(z, jjTf) is set to be the vacuum fragmentation func- 
tion whereas da me ^^ x /dP T dy is computed from the 
same equation with Df^h(z, pA) replace by the suitably 

averaged MMFF {Dmm{z, M/))t.aa' This averaging has 
to be done over all possible paths of partons through the 
medium. 

The probability density P(xo,yo) for finding a hard ver- 
tex at the transverse position ro = (xo,yo) and impact 
parameter b is, again in leading order, given by the prod- 

I 



Using the approximate scaling relation described in sec- 
tion [TII3 the medium modified fragmentation function 
Dmm{z, /i 2 , C) f° r a path £ can be found by computing 
the line integrals / d(q(( or f d(D(() over Eqs. (fl7 ) .(fig] ) 
or by evaluating Eq. ([2H)l respectively. The MC shower 
code is then used to compute Dmm(z, p 2 ,C) f° r eacn 
value of AQ 2 ot , AE tot or f med obtained. 
As discussed in more detail in ,19(, there is a concep- 
tual problem with using a MMFF computed for a fixed 
partonic scale in Eq. (f2"Tj) where D(z : ^) is an object 
defined at a given hadronic scale. This is a generic prob- 
lem in obtaining fragmentations from a MC code which 
starts with given parton properties, however in practice 
the scale evolution in the RHIC kinematic range is small 
as compared to other uncertainties in the computation 
and the resulting uncertainty can be tolerated. In the 
following, we use a MMFF determined at the partonic 
scale \i = 20 GeV. 



B. Comparison with data 

With the medium given by the hydrodynamical evolution 
model described in [34| and the expressions for hadron 
production in vacuum and medium Eq. (|21[) . the remain- 
ing unknown quantities for a comparison with data are 
the parameters K, Kr> and Kf which link the medium 
properties in terms of the energy density e with the 
parton- medium interaction parameters q, D and fmed- 
Note that according to the results of IIVDI the value of 
these parameters cannot be uniquely determined for sin- 
gle inclusive hadron spectra, but depends on the choice 
of the scale Q m in in the shower simulation. In the follow- 
ing, we show the best fit of K, Kr> and Kf to the data 
given the choice Q m in = 1 GeV. 

In Fig. [5] we show the calculated nuclear suppression fac- 
tor as a function of hadron momentum Pt for all three 
scenarios in comparison with the data for 7r° production 



uct of the nuclear profile functions as 

p , x T^(r +b/2)r A (r -b/2) 

P(x , y ) = , (27) 

where the thickness function is given in terms of 
Woods-Saxon the nuclear density pa(t,z) as Ta(i-) = 
/ dzp A {r, z). The MMFF must then be averaged over 
this quantity and all possible directions <f> partons could 
travel from a vertex as 



(28) 
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FIG. 9: Calculated nuclear suppression factor for three differ- 
ent scenarios for the parton-medium interaction (see text) as 
a function of hadron momentum Pt in comparison with data 
by the PHENIX collaboration |7|. In all cases, the relevant 
constant relating medium energy density and parton-medium 
interaction parameter has been fit to data. 



in 200 AGeV central Au-Au collisions. The most strik- 
ing observation is that all three scenarios are surprisingly 
similar and could not possibly be distinguished by cur- 
rent data for Raa- Most notably, all scenarios exhibit a 
falling trend as Pt where scenarios based on leading par- 
ton energy loss typically exhibit a rising trend (see e.g. 
[HIH). In [H, this property was tentatively attributed 
to the fact that a description of the whole shower keeps 
track of multiple soft hadron production. It was also sug- 
gested that the same physics underlies the enhancement 
of dN/ d£ in the large £ region and the falling of Raa with 
Pt- However, the present investigation shows that both 
ideas must be discarded, as the DRAG scenario in which 
no enhanced soft hadron production occurs shows also no 
enhancement of dN/d£, but the same falling trend of Raa 
with Pt as the other scenarios. Thus, the falling trend is 



7/>27T /• OO /• OO 





(Dmm(z,h 2 ))t aa =-^- I d<p I dx I dy P(xo,yo)DMM(z,n 2 ,0- 
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not a phenomenon characteristic of radiative energy loss 
but substantially more general. It also has been observed 
in other models where a modification of the whole shower 
by the medium was considered, cf. e.g. (25l.l49j. 



duced into the simulation of in-medium shower evolution. 



C. Extracting medium parameters 
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FIG. 10: Comparison of the MMFF of a 20 GeV d-quark into 
charged hadrons for a schematic ASW leading parton energy 
loss scenario and two different scenarios in which the whole 
shower is evolves (see text). All parameters are chosen such 
that the curves agree at z — 0.6. 

In order to gain greater insight into the differences be- 
tween scenarios which compute energy loss for the leading 
parton and between those where the whole shower evolu- 
tion is modified by the medium, we present a schematic 
comparison between the MMFFs in the RAD, the DRAG 
and the ASW scenario in Fig. [TU] (the FMED scenario, 
being in essence indistinguishable from the RAD scenario 
is not shown here). 

It has been pointed out repeatedly (see e.g. [13, EH) 
that the ASW scenario applied to RHIC kinematics in 
essence leads to complete absorption of about 75% of all 
partons, ~ 15% emerge without any energy loss and only 
a small fraction is found after finite energy loss. To good 
approximation, the MMFF in the ASW scenario is thus 
just a downward shift of the vacuum baseline. 
This has been done in Fig. 1101 where all parameters have 
been adjusted such that the curves agree at z — 0.6. It 
is obvious that the shape of the schematic ASW result 
is quite different fron the other scenarios. In particular, 
the difference between RAD and DRAG is much less pro- 
nounced than between either of those and ASW. It is in 
essence given by the presence or absence of soft hadron 
production and confined to the region z < 0.2. Thus, it 
appears that the different curvature of D(z) at z > 0.5 
is responsible for the rising vs. falling trend in Raa > and 
thus the way the high Pt end of the shower evolves rather 
than low Pt hadron production are seen in the data. 
At present, the falling trend seems not to be supported 
by the data. Should this be confirmed by more precise 
measurements, presumably the possibility of complete 
absorptions of partons by the medium needs to be intro- 



We can use the above results to tentatively extract 
medium properties. In the RAD scenario, K is a di- 
mensionless parameter and from the fit shown in Fig. [5] 
the value K = 3 is found. This differs from the result in 
[l|| where K — 1.5 was obtained, note however that in 
the present work a randomziation of the formation time 
(see EQ. (fTTjl ) has been performed and that the effect of 
the medium prior to thermalization has been included in 
a schematic way. These two differences account for the 
changed value of K. 

With this value of K, qo, i.e. the highest transport co- 
efficient reached in the evolution in the medium cen- 
ter at thermalization time of 0.6 fm/c is found to be 
15.6 GeV 2 /fm when Q m i n — 1 GeV is assumed. For 
Qmin — 0.7 GeV, the extracted value olK changes to 1.4 
and q = 7.2 GeV 2 /fm. 

Kn is a dimensionful parameter which can be expressed 
in units GeV -1 . The same fit yields (due to the numerical 
coincidence mentioned before) Dq = 15.6 GeV/fm for 
Qmin = 1 GeV and D = 7.2 GeV/fm for Q mm = 0.7 
GeV. 

While these numbers appear large, it has to be remem- 
bered that they reflect a snapshot of the medium at its 
peak density, from which the energy density drops rapidly 
as a function of time due to the expansion. Since f me d is 
not in connected to any microscopical properties of the 
medium, we refrain from analyzing its value here. 



VI. DISCUSSION 

We have presented a comparison study of three differ- 
ent mechanisms for the parton-medium interaction in the 
framework of an in-medium shower evolution. In this 
study, we have considered a variety of assumptions about 
the evolution of the medium, among them a constant 
medium with different length L, an evolving medium for 
parths from the medium center, paths from the medium 
surface and an average over all possible paths in the 
medium. We have considered three different types of 
shower initiators — heavy quarks, light quarks and glu- 
ons. We have studied single parton observables such as 
the distribution of the leading parton momentum or the 
energy loss probability density P(AE), single hadron ob- 
servables like Raa as well as multihadron observables 
such as the hump-backed plateau dN/d£. In addition, 
we have also studied merely technical aspects of mod- 
elling such as the role of the cutoff parameter Q m in or 
the effect of randomizing the formation time of partons 
in branching. From the results in all these different situ- 
ations, some generic properties can be identified. 
• The nuclear suppression factor Raa is not a good ob- 
servable to distinguish different scenarios of the micro- 
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scopical interaction of partons with the medium. This 
statement has been made previously from different angles 
(see e.g. [H, [H3,[5H) and the present results merely con- 
firm previous findings in yet another framework. More 
differential observables are needed to determine the na- 
ture of parton-medium interaction. 

• The falling trend of Raa as a function of Pt is appar- 
ently unrelated to low Pt multi hadron production and 
rather a generic feature observed in models which do not 
consider energy loss from a leading parton but rather a 
modification of the whole shower. For example, the re- 
sults of the Higher Twist approach applied to the leading 
parton show a rising trend of Raa with Pt [52| , however 
when resummed in the shower evolution equations and 
applied to the whole shower, the Higher Twist approach 
leads to a falling trend [3(|, i.e. the same observation 
is made in quite a different framework. If future data 
confirm a rising trend, non-trivial modifications to the 
shower evolution codes, such as the possibility of com- 
plete parton absorption by the medium, need to be con- 
sidered. 

• The properties of medium-induced radiation as a mech- 
anism for the parton-medium interaction appear rather 
generic. There is no observable in this study in which 
the RAD and the FMED scenarios lead to substantially 
different results. The useful implication would be that 
in many observables it is really the underlying physics 
mechanism one is probing, not technical details of how 
this mechanism is implemented in a particular model. 

• In contrast, a different physics mechanism as exem- 
plified here by the DRAG scenario appears distinct in 
several quantities. Not only is its excitation function in 
terms of mean energy loss as a function of medium den- 
sity different than for radiative scenarios (which could be 
tested by variations in collision centrality), but also the 
absence of soft hadron production induces pronounced ef- 
fects in jet observables such as the angular distribution of 
hadrons around the jet axis or the hump-backed plateau. 
However, the need to identify a jet in a heavy- ion collision 
above the soft background introduces additional compli- 
cations. In essence, a medium-modified jet has properties 
different from a jet in vacuum and is hence less likely to 
be identified as jet. A measurement of jets must be care- 
fully designed to avoid this trigger bias which tends to 
hide the very effect one would like to study [2(| ■ 

• Technical aspects of the modelling, such as the choice 
of Qmin or the randomization of the formation time as 
investigated here, do not appear to change the results 
qualitatively. However, there is a substantial ambiguity 
once one tries to extract quantitative medium parameters 
from the computation, especially when this extaction is 
based on a single observable. As is the case for vacuum 
shower codes, the relevant technical model parameters 
should eventually be determined by the best fit to a large 
body of data. 

There are several more properties of the parton-medium 
interaction which could be exploited to distinguish dif- 



ferent scenarios. A very promising candidate is the path- 
length dependence of the medium effect. Experimentally, 
this can be varied moderately by considering the nuclear 
suppression as a function of the angle of hard hadron 
with the reaction plane [lOj or more strongly by consid- 
ering back-to-back correlations (which however require 
a careful modelling, as the relevant geometry arises as a 
complicated function of the geometrical bias of the energy 
loss on the trigger hadron itself [H, Q3] • For example, in 
(33| is was argued based on the different pathlength de- 
pendence that elastic energy loss cannot be responsible 
for the suppression of hard back-to-back hadron correla- 
tions. 

In the present paper, we have refrained from making any 
comparison based on pathlength dependence. Such a 
study (which is quite substantial on its own) along with 
a comparison with the experimental results on back-to- 
back correlations and the variation of the suppression as 
a function of the reaction plane angle will be the topic of 
a future publication. 



VII. OUTLOOK 

The study presented here shows that jet observables are 
more powerful in order to distinguish different microscop- 
ical physics process of the parton-medium interaction 
than observables which are only sensitive to the lead- 
ing hadron. However, the need to identify a jet above 
the background medium may quickly eliminate this ad- 
vantage, at which point one has to resort to what has 
been termed the 'golden channel' — 7-jet correlations in 
which the presence of a photon not only allows to get an 
unbiased jet sample but also reveals the kinematic of the 
jet. Unfortunately, due to the smallness of the electro- 
magnetic coupling, the statistics in this channel is poor 
and the measurement is difficult. 

Future jet measurements at RHIC may overcome this 
problem by high luminosity, whereas future measure- 
ments at LHC where the scale separation between a hard 
process and the soft medium is considerably larger than 
at RHIC may not suffer significantly from trigger bias at 
all. 

There is now good reason to assume that jet observables 
will reveal important information about the microscopi- 
cal properties of the medium, and Monte Carlo simula- 
tions of in-medium showers such as YaJEM or JEWEL 
will most likely be the appropriate tools to extract this 
information. 
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